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Abstract 

The production of charged Higgs boson pairs from gluon- 
gluon fusion is studied in the minimal supersymmetric 
model(MSSM) at proton-proton colliders. We find that the 
contribution of the scalar quark loops to the subprocess pro- 
duction rate is substantial and consequently the production 
rate in proton-proton colliders is enhanced in the MSSM if 
there exist scalar quarks with favorable masses. The results 
show that this cross section may reach a few femptobarn in 
the future LHC collider with plausible values of the param- 
eters. 
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1. Introduction 



The minimal standard model(MSM) |[T|[0], a gauge theory based on SU{3)c x SU{2)l x 
U{1)y, is a very successful model which nearly describes either quantitatively or qualitatively 
all available data pertaining to the strong and electroweak interaction phenomena. Until now 
the symmetric breaking sector of the electroweak interactions has not yet been directly tested 
experimentally. The future multi-TeV hadron colliders such as the CERN Large Hadron 
Collider (LHC) and Next Linear Collider (NLC) are elaborately designed in order to detect 
a 'light' scalar boson(mj:/ < ITeV) associated with the symmetry-breaking mechanism. 

Any enlargement of the sector beyond the single SU{2)l Higgs doublet of MSM, with 
two or more doublets as required in supersymmetric theory, necessarily involves new physical 
particles. In this extension theory the physical spectrum includes charged Higgs bosons. 
Like the general two-Higgs-doublet model(THDM), the minimal supersymmetric standard 
model(MSSM) also requires two SU{2) doublets to give masses to up and down quarks |@|@]. 
Since both models contain the same Higgs sector structure, it is of interest to determine 
whether THDM or MSSM Higgs sector might appear in future experiments. 

With regard to colliders in the future there are several mechanisms which can produce 
charged Higgs bosons. Among them the charged Higgs boson pair productions via gluon- 
gluon collision, quark-antiquark annihilation and 77 collision, play a significant role in study- 
ing charged Higgs bosons. Similar to the neutral Higgs pair production process in hadron 
colliders, the charged Higgs pair production process also involves Higgs sector interactions, 
therefore the measurement of these kinds of couplings could be a good test of our theory. 
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Recently Plehn et al provided calculations concerning the pair production of neutral Higgs 
particles in gluon-gluon collisions in MSSM[|] and charged Higgs boson pair production at 
NLC in 77 modes are also investigated in references 0. The charged Higgs boson pair 
production via gluon-gluon fusion in general THDM is also represented in Ref.0. In this 
reference, it is concluded that the production rate of this process in THDM is not sensitive 
to the neutral Higgs boson masses and the couplings of self interactions of Higgs bosons, but 
it is strongly related to the couplings of the charged Higgs boson to quarks. Willenbrock 
discussed the charged Higgs boson pair production via gluon fusion in MSSM in Ref.0. He 
showed that the gluon fusion via quark and scalar-quark loops would dominate the usual 
electroweak quark pair annihilation process if there exist sufficiently heavy quarks or scalar- 
quarks. In his work, he took some approximations, such as the large quark and squark mass 
limits, i.e., mg(m|) >> s, in numerical calculations. After we have finished our paper, we 
found a work also on this process. One can refer to Ref. ||1CI||. 



In this paper we concentrate on the capability of the charged Higgs boson pair production 
via gluon-gluon collisions at pp colliders in frame of MSSM without the approximations used 
in reference 0. The paper is organized as follows: The analytical evaluations are given in 
section 2. In section 3 there are numerical results, discussion and a short summary. Finally 
the explicit expressions of the relevant form factors induced by the supersymmetric quark 
loop diagrams are collected in Appendix. 

2. Calculation 
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The MSSM theory gives a large number of parameters, since it predicts more than dou- 
bhng the MSM spectrum of particle states. The number of the parameters can be reduced 
largely by embedding the low-energy supersymmetric theory into a grand unified(GUT) 
framework. Besides the parameters which can be achieved in supergravity model[Q, we need 
extra parameters to describe the Higgs sector: one Higgs mass parameter and the two vac- 
uum expectation values of the neutral Higgs fields. We can also use tan/3 = ^, the ratio of 
the vacuum expectation values of the two-Higgs-doublet fields which break the electroweak 
symmetry. If we ignore the flavor mixing, the mass matrix of a scalar quark takes the 



following form||ll||: 



-^-('1 'ocl mirj w 



here and are the current eigenstates and for the up-type scalar quarks. 



=Ml+ml cos 2/5(^ - sin Ow^) + mj. 



m?^ = MIj + m| cos 2(3Qq sin Ow"^ + m^q-, 

aq = fi cot (3 + A^M, (2) 

for the down-type scalar quarks, 

m?^ = - m| cos 2/?(^ + Qq sin Ow^) + mj, 

m?^ = + m\ cos 2/3(5g sin Ow"^ + w-g, 

aq = iitixn(3 + A*M, (3) 
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where Qq is the charge of the scalar quark, Mq, and are the self-supersymmetry- 
breaking mass terms for the left-handed and right-handed scalar quarks. x M is a trilinear 
scalar interaction parameter, and is the supersymmetric mass mixing term of the Higgs 
boson. When there is no CP violation, we can regard that is a real constant. At the 
Plank mass scale we take 

Ml^Ml^Ml^ (4) 
The mass eigenstates q\ and q2 are related to the current eigenstates qi and qu by 

qi = qt cos 9q + qn sin 9q 

h = -qL sin 9q + qn cos 9q. (5) 

and 



2aqmq 

tan26'5 = ^ — (6) 



The masses of qi and g2 are 



In the evaluation we neglected the first SUSY generation mixing angle, whereas in the 
second SUSY generation we took sin^^ = 0.2, and let the third SUSY generation has the 
mixing with sin^^ k, due to the large splitting of the mass eigenstates. The numerical 
value of the lightest physical stop quark mass is taken as m^^ — 180GeV and the sbottom 



quark mass as m^^ = 4006*6^. We set the masses of the squarks of ui, di, ci and si 
degenerated with mass value SOOGeV^. As a quantitative example, we took /i = in our 
calculation. 

The generic Feynman diagrams contributing to the subprocess gg H^H^ which are 
involved in frame of MSSM are depicted in figure 1. The diagrams created by exchanging 
the two external gluon lines are not shown in Fig.l. Fig.l(l ~ 10) are the same as those in 
the THDM model 0, we call them as "THDM part". Fig. 1(11 ~ 32) are the diagrams with 
squark-loop. There is no tree level contribution for this subprocess, therefore the proper 
vertex counterterm cancels with the counterterms of the external legs diagrams. That is to 
say that the evolution can be simply carry out by summing all unrenormalized reducible and 
irreducible diagrams. The result is finite and gauge invariant. The relevant Feynman rules 



can be found in references |TT|] [jT2|. Since the main contributions are due to the loops of heavy 
quarks and the scalar quarks, we neglect the light quark loop diagrams. We can see from 
Fig.l that there are mainly three kinds of production mechanisms. One is by exchanging 
a photon or a boson to produce a pair of charged Higgs bosons, we refer to this as 
7 - Z° exchange s-channel (see Fig.l(l,2), Fig.l(ll,12), Fig.l(15,16) and Fig.l(19,20)). For 
the 7 — Z° exchange s-channel diagrams, the contributions are shown to be zero. This is 
the consequence of Furry theorem. Since the Furry theorem forbids the production of the 
spin-one components of the and 7, and the spin-zero component of the Z^ vector boson 
does not couple to a pair of scalar particles with the same mass. The second mechanism 
is through producing the virtual neutral Higgs and bosons which are coupled to 
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gluons by quartic coupling including squarks or by a triangle of heavy quarks or squarks, 
then the charged Higgs boson pair is appeared by the subsequent decay of the virtual Higgs 
bosons (see Fig.l(3,4,13,14,17,18)). It is called neutral Higgs boson exchange s-channel. An 
alternative is that the charged Higgs bosons are produced by means of the box and quartic 
interactions. 

In this paper we follow the notations in reference that pi, p2, ps and represent 
the four-momenta of the final charged Higgs bosons and the initial gluons, respectively. We 
write the corresponding matrix element for each of the diagrams in Fig.l in the form as 

Ml = Y^3s^ti{P3)^u{p4)A2 ■ 

ihig'" + f2,M + + kiv'2v\ + hip'2P2 + 

/ll,^e''"%i,aP2,/3P3,7P^ + /l2,/e'^"%i,aP2,/3P3,7P2)- (§) 

We denote Afc(k=l,2) as the color factors of the diagrams with and without gluon-gluon- 
squark-squark quartic vertex respectively and it can be expressed as 



Al = Tr[^5aJ + dabcT'']=5ab. 

A2 = Tr[T^T'] = ^6ab. (9) 

In above equations / is a 3 x 3 unit matrix. Qs is the strong coupling constant and 
T"(a = 1 ~ 8) are the SU{3)c generators introduced by Gell-Mann and dabc are the structure 



const ant s||TB||. Here we should declare that in Eq.(8) A2 = jSab is the common factor of the 
all color factors in amplitudes and the remained parts of color factors are collected in form 
factors. 

The total matrix element of subprocess gg H^H^ in MSSM can be expressed as 

10 32 

1=1 n=ll l-LUJ 

= Mthdm + MsuSY, 

where Mthdm represents the summation of the amplitudes of the THDM part and Msusy 
is the total amplitudes of squark-loop diagrams. 

Then we can split each form factor appeared in total matrix element up into two parts. 

/. = /f™^ + /r"^(^ = l~12). (11) 

32 10 

where ff^''^^{i = 1 ~ 12) = E and ff"^^'{i = 1 ~ 12) = E fi,i- The explicit expres- 

n=ll ' 1=1 

sions of the form factors fi^^^{i = 1 ~ 12), which come from the squark-loop diagrams, 
are listed in Appendix. For the form factors of the THDM part (i.e. f^^^'^\i = 1 ~ 12)) 
one can find them in reference|0]. The total cross section for subprocess gg H^H~ can 
finally be written in the form 

where \M\ is the initial spin and color averaged matrix element squared and = {m'j^± — 
^s) ± ^sPh±- The total cross section for the charged Higgs pair production through gluon- 
gluon fusion in proton-proton collisions can be obtained by integrating the a over the gluon 
luminosity. 

a{pp ^gg^ H+R- + X) = C dr^^a^gg H+H~ at s = ts). (13) 

J4m2 , /s (It 



where a/s and denote the proton-proton and gluon-gluon cm. energies respectively 



and ^^^4^ is the gluon luminosity, which is defined as 



dCgg dxi 



(14) 



Here we used r = X1X2, and one can find the definitions of Xi and X2 in Ref.|0]. In our 
calculation we adopt the MRS set G parton distribution function Fg{x) [|T1|, and ignore the 
supersymmetric QCD corrections to the parton distribution functions for simplicity. The 
numerical calculation is carried out for the LHC at the energy around 10 ~ lATeV. 

3. Numerical results and discussion 



For the numerical evaluation we take the input parameters ||15|| as rrii, = 4.5 GeV, mz 



91.1887 GeV, mw = 80.2226 GeV,mt = 175 GeV , Gp = 1.166392 x 10-^{GeV)-^ and a = 
iffo^- We adopt a simple one- loop formula for the running strong coupling constant as as 

««(^) = , , 33-2n, , / ^ V (15) 

where as{mz) = 0.117 and Uf is the number of active fiavors at energy scale /i. The two 
representative quantities of the two vacuum expectation values tan/3 =1.5 and tan/? = 30 
are chosen. 

We take the mass of charged Higgs bosons mH±, and the ratio of vacuum expectation 
values as free parameters in MSSM. The other Higgs boson masses and the mixing angle a 
are fixed by quoting the MSSM relations in terms of tan/? and mH±- The relevant formulas 
are given as 
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„ TTi A H~ Tfi '7 2 2 2 //^\ 

tanza = tanzp — 5 5-, where 'in a = m[j± — mw. (16) 

— rriz 



2 _ 1 

'^//o./iO — 2 



(17) 



The cross sections of the subprocess gg — > H^H" with mij± = 150 GeV are depicted 
in Fig.2(l), on the conditions of the mass vahies of squarks as mentioned above. In this 
figure, the small peak around a/I ~ 350 GeV of each curve originates from the resonance 
effect, where -\/I ~ 2mt and ~ 2771^^ are satisfied. Around -\/I ~ 800 GeV , there is 
another higher peak for each curve, which comes from the enhancement of the resonance 
effects where the cm. energy of gluon-gluon system is located at -\/I ~ 2mf2 ~ 800 GeV and 

~ 2m^^ ~ 800 GeV . Fig. 2 (2) shows the cross section of this subprocess as the functions 
of cm. energy of incoming gluons with the same neutral Higgs boson, quark and squark 
masses as in Fig. 2(1), but mH± — 300 GeV. In Fig. 2(2) again the resonance effects of the 
squarks are obvious at the position near \/I - 800 GeV. 

The dependence of the subprocess cross section to the ratio of vacuum expectation values, 
is plotted in Fig. 3 with mH+ = ISOGeF. The two curves correspond to a/I = AQQGeV and 
-\/I = 800Gey , respectively. They all have the minimal cross sections when tan {3 is around 
6.5. The figure shows that the cross section of the subprocess increases when the a/I goes 
up from 400 GeV to 800 GeV . The dependence of the subprocess cross section to the mass 
of charged Higgs bosons is shown in Fig.4 with a/I = ITeV. The cross section increases 
in the region where the charged Higgs boson mass has the value approximately larger than 
~ 300 GeV. That is because when the value of the the charged Higgs boson mass approaches 
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|\/^ = 500 GeV, the production rate will be enhanced by the threshold effect, but when the 
mass of charged Higgs boson is very close to 500GeV, the cross section drops down quickly 
due to the phase space suppression. The cross sections of process pp ^ gg ^ H'^H" + X in 
the LHC energy regions in both MSSM and THDM are depicted in Fig.5(l~2). The charged 
Higgs boson pair production rates in MSSM are read to be about 1.9 to 5.2 femptobarn at 
the LHC energy regions when mH+ = l50GeV. With the heavier charged Higgs mass, the 
production rate will be smaller. The production rates in the MSSM are higher than those 
in the THDM with the same parameters. In our calculation, we also find that among the 
squark-loop diagrams, the third generation squarks make the decisive contribution, while 
the contribution from the other squarks is rather small comparatively due to the decoupling 
of the heavy squarks. 

In conclusion, we investigated the pair production process of charged Higgs bosons via 
gluon-gluon fusion in pp collider at LHC. The numerical analysis of its cross section is carried 
out in the MSSM model. With the possible parameters, the cross sections at future LHC 
collider may reach 5.2 femptobarn. The calculation shows that the contribution from the 
squark-loop diagrams is obvious and consequently enhances the cross section of pp ^ gg ^ 
H^H~ + X, if there exist scalar quarks with favorable masses. 

We express our thanks to Dr. M. Spira for useful discussions. This work was supported 
by the National Natural Science Foundation of China. Part of this work was done when 
two of the authors. Ma Wen-Gan and Yu Zeng-Hui, visited the University Vienna under an 
exchange agreement (project number: IV.B.12). 
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Appendix 

The form factors involved in equation (8) are presented explicitly as: 

fl = -2^5 XI { Yl {ghOH^H^ghOqiqAhO + gHOH±H^gHOqiqAHO + gH+H-Mi) 

{4C24[-P3,P1 + P2, mq.,mq^,mq.] - Bo[-pi - P2 , 777,^,, 771^,]} 

{ci,Sj),{ti,bj) 

J=l,2 (qy.^i,qdj)=(ui,dj) 
- ^C'obl + P2, -P2, mq^,vmq^,vmq^,^] 
+ iD27 [-P3 ,Pl+P2, -P2 , rUq^ , , rUq^ . , TJlq^ . , TUq^ . ] 

+ iD27 [-P3 , Pi + P2 , -P2 , rrig^j , rrig^j , rn^-^^ . , nig^ .] 

+ iD27 [pa , -Pi - P2 , Pi , rrig^ ^ , rrig^ . , rUg^ . , Tn,^ . ] 

+ ^-027 , -Pi - P2 , Pi , rrig^^i , rug^ . , rriq^ . , .] 

+ D27 bi , P2 - Pa , -P2 , rrig^j , rrig^ . , r77g-,_^^ , m^-^ J 

+ ^27bi>P2 - Pa, -P2, mg^,i,mq^j,'rnq^,i, rriq^J}} - 1) 



f2 = -Sfifg X { XI {ghOH^H^fghOMi'^hO + gHOH^H^FgHOqiqiAffO + QH+R-qiqi) 

(C23 - (^22) [-Pa, Pi + P2, mq^,mq^, rriq^] 

(ci,5j),(ij,6j) 

+ E E _ ^l.±,.,,,,,{2i(^24 

-D22)[-P3,P1 +P2, -P2, 
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+ 2i{D24 - D22) [-P3, Pi + P2, -P2, . , ruq^ . , mg-^_ . , m,^ 

+ 2i(D24 + 2D26 - £>22 - -D23 - -D25)[P3, "Pi " P2,Pl, , , , 

+ 2^(D24 + 2^26 - £'22 - -D23 - D25) bs, -Pi - P2, Pi, m,^,,- , m^dj ' ™?dj ' ™««,J 

- (i:>o + 3Dii + 2i:>2i)[pi,p2 - P3, -P2, mq^^pmq^,i,mq^,prnq^^ 

- {Do + 3Du + 2D2i)\pi,P2 - P3, -P2, mg^,i,mg^j,'rng^,i, mg^J}} - 2) 



fs = -25^^ J2{ J2 {ghOH±H^gh°qiqAhO + QH" H±H^ gHOqiqiAnO + gH+H-Mi) 

(C23 - C22) [-P3, Pi + P2, rug^ , mg, , nig^] 

+ E E _ 9%^q.,q,AnD24 + D2e 

i—^'"^ (qu,i,qd,j)={ui,dj) 

-^22)[-P3,Pi + P2, -P2, mg^,i,mg^,i,mg^,i,mg^J 
+ 2i{D24 + 1^26 - ^22) [-P3, Pi + P2, -P2, nig^^^ , uig^^^ , nig^^^ , m^-^ J 

+ 2i(D24 + D2& - D22 - -D25)b3> -Pl - P2,Pl, mg^,vT^qu,vT^^u,V^lid,j\ 

+ 2i ( D24 + -D26 - ^22 - ^25 ) , -Pi - P2 , Pl , . , nig^ . , nig^ . , rrig^ .] 

- (Do + Dn + 2Di2 - 2D13 + 2D24 - 2D25)bi>P2 - Ps, -P2, mq^j,mq^,i,mq^,j,mq^,i] 

- (Do + Dn + 2Di2 - 2D13 + 2D24 - 2D25){pi,P2 - Ps, -P2, mq^,i,mq^j,mq^,i,mq^J}} {A - 3) 



/4 = -2fifs E { E {ghOH±H^ghOqiqi^hO + gHOH^H^gHOqiqiAnO + gn+H-qiqi) 

^=^>'^ qi=Ui,di,Ci 

(C23 - C22) [-P3, Pl + P2, ruq. , m^, , rriq^] 
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-D22)[-P3,P1 + P2, -P2, mq^,i,mq^,i,mq^,i,mq^J 
+ 2i{D24 + D26 - D22)[-P3,P1 + P2, -P2 , niq^. , Ulq^. , TUq^. , , ] 

+ 2i(D24 + D2Q - D22 - -D25)[P3, -Pl -P2,Pl, mq^,i,mq^,i,mq^,i,mq^J 

+ 2i{D24 + D26- D22 - D25 ) ba , -Pi -P2,Pi, ruq^ . , rUq^ . , ruq^ . ,mq^ .] 

- (D12 - Dis + 2D24 - 2D25)\pi,P2 - P3, -P2, mqjj,mq^,i,mq^j,mq^,i] 

- (D12 - Di3 + 2D24 - 2D25)Ipi,P2 - P3, -P2, mq^,i,mg^,j,mg^,i,'mq^J}} - 4) 



h = -2fi's E { E {9h0H±Hl^gh0qiqAh0 + 9H0H±H1^9H0qiqiAH0 + gH+H-qiqi) 

(C23 - C22)[-P3,P1 + P2, mq^,mq^,mq^] 
{ci,Sj),{ii,bj) 

+ E E , ^1-,-...-.. {2^(^24 

■'"-^'^ {qu,i,qd,j)={ui,dj) 

-D22)[-P3,pi + P2, -P2, mg^,i,'mq^,i,mg^,i, rrig^J 

+ 2i{D24 - D22)[-P3,P1 + P2, -P2, rJlg^. , TTlq^. , TJlg^. , TTlq^J 

+ 2i(D24 + 2^26 - D22 - D23 - -D25)[P3, "Pl - P2,Pi,mq^ ^,mq^ .,mq^ .,mq^ .] 

+ 2i{D24 + 2D2Q - D22 - D23 - D25)\P3, -Pl - P2,Pl, mq^,j,mq^j,mq^j,mq^,i] 

- (D12 - Di3 + 2D22 + 2D23 - 4.D2e)\pi,p2 - P3, -P2, mq^,j,mq^,i,mq^j,mq^,i] 

- (D12 - Di3 + 2D22 + 2^23 - 4:D2e)\pi,P2 - P3, -P2, mq^,i,mq^j,mq^,i,mq^J}} " 5) 
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/6 = /7 = ■ ■ ■ = /l2 = 



(A -6) 



= 2 — r~T^ ' 

s — m^o + li HorriHO 



AhO = . {A -7) 



The coupling constants are obtained by converting Feynman rules for Xgg(where X is a one- 



or two-particle state) from the — qji basis, their Feynman rules can be found in Ref . 
to the qi — q2 basis. For the converting formulae one can refer to reference[Q. 

In above expressions we adopted the definitions of one-loop integral functions in reference 
and defined d = 4 — 2e. The arguments of two-point, three-point and four-point integral 
functions are written at the end of formulae in brackets. The numerical calculation of the 
vector and tensor loop integral functions can be traced back to four scalar loop integrals Aq, 



Bq, Co, Dq as shown in the reference 117 
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Fig.l The Feynman diagrams of subprocess gg — > H^H~ . Figures (1 ~ 10) are so called 
THDM diagrams. Figures (11 ~ 32) are the diagrams having one squark-loop. (1 ~ 4, 
11 ~ 20, 27 ~ 32) are s-channel diagrams. The others are box diagrams which include t- 
and u-channels. 

Fig. 2(1) Total cross sections of the subprocess gg — > H^H~ as the functions of the 
with mH± = 150 GeV in the MSSM. The full-line is for tan/? = 1.5. The dotted-line is for 
tan/? = 30 . 

Fig.2(2) Total cross sections of the subprocess gg — > H^H~ as the functions of the 
with mH± = 300 GeV in the MSSM. The full-line is for tan/3 = 1.5. The dotted-line is for 
tan/? = 30. 

Fig.3 Total cross sections of the subprocess gg — > H^H~ as the functions of the tan/? with 
mH± = 150 GeV in the MSSM. The full-line and dotted-hne correspond to = 400 GeV 
and \/S = 800 GeV respectively. 

Fig. 4 Total cross sections of the subprocess gg — > H^H~ as the functions of the mH± with 
= 1 TeV in MSSM. The fuU-hne is for tan/3 = 1.5. The dotted-hne is for tan/3 = 30. 

Fig. 5(1) Total cross sections of the process pp ^ gg ^ H^H~ -|- X as the functions of the 
i/i with mH± — 150 GeV at the LHC energy region. The fuU-hne is for tan/3 = 1.5 in the 
MSSM. The dashed- line is for tan/3 = 1.5 in the THDM. The dotted-line is for tan/3 = 30 
in the MSSM. The dashed-dotted-hne is for tan /3 = 30 in the THDM. 



18 



Fig. 5 (2) Total cross sections of the process pp ^ gg ^ H^H + X as tlie functions of the 
with mH± — 300 GeV at the LHC energy region. The fuU-hne is for tan/3 = 1.5 in the 
MSSM. The dashed-hne is for tan /? = 1.5 in the THDM. The dotted-hne is for tan /3 = 30 
in the MSSM. The dashed-dotted-hne is for tan/? = 30 in the THDM. 
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